Nonlinear statistics of quantum transport in chaotic cavities 
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In the framework of the random matrix approach, we apply the theory of Selberg's integral to problems of 
quantum transport in chaotic cavities. All the moments of transmission eigenvalues are calculated analytically 
up to the fourth order. As a result, we derive exact explicit expressions for the skewness and kurtosis of the 
conductance and transmitted charge as well as for the variance of the shot-noise power in chaotic cavities. The 
obtained results are generally valid at arbitrary numbers of propagating channels in the two attached leads. In 
the particular limit of large (and equal) channel numbers, the shot-noise variance attends the universal value 
1/64/3 that determines a universal Gaussian statistics of shot-noise fluctuations in this case. 
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Quantum transport of non-interacting electrons in meso- 
scopic systems can be conventionally described in the frame- 
work of scattering theory. 1,2 The conductance and shot-noise 
are the brightest and, perhaps, most frequently considered ex- 
amples. Being expressed in terms of transmission eigenvalues 
Ti of a conductor, the dimensionless conductance, g, and the 
zero-frequency shot-noise power, p, of the two-terminal setup 
at zero temperature are respectively given by 



^T, and p = Y J T l {l-T l ). 



(1) 



Here, n = min(/Vi, N%) where is the number of scatter- 
ing channels in each of the two attached leads. 

Generally, the positive quantities < 1 are mutually cor- 
related random numbers whose fluctuations depend on the 
conductor's nature. In the case of chaotic cavities considered 
below, the joint distribution of T s ; is believed to be provided 
by the random matrix theory (RMT) and reads as follows 1 



mm = Afp 1 n t?- 1 n i T i - T ^ ■ 



(2) 



j<k 



with a = ^(\Ni — 7V2I + I) and normalization constant Mp. 
The symmetry index (3 (= 1, 2, or 4) distinguishes between 
the three standard RMT classes (orthogonal, unitary, or sym- 
plectic ensembles, respectively) which are realized depending 
on the presence or absence of time-reversal and spin-flip sym- 
metry in the system.^ 3 - 

The exact RMT results for the average and variance of the 
conductance are well known for quite a long time 1,4,5 and that 
for the average shot-noise power has become available only 
recently^ As concerns higher order cumulants of these quan- 
tities, their exact RMT expressions valid at arbitrary N\ $ have 
not been reported in the literature so far, a progress in this di- 
rection being announced very recently^ 

In this communication, we answer this question by further 
developing and applying the theory of Selberg's integral to 
the problem. As was recently recognized, 6 such an approach 
is a powerful non-perturbative method suited particularly well 
for the studies of moments and counting statistics, see also 
Refs. 7 and 8. It represents a useful alternative to orthogonal 



polynomial 3 or diagrammatic 9 approaches, especially in the 
situation when the finite dimensionality of relevant random 
matrices becomes important. 

Selberg's integral appears naturally in the problem first as 
an integral determining the normalization constant 
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This expression assures that (O is a probability density, being 
generally valid for discrete n and continuous a and To 
study the cumulants of g and p, one needs to know the mo- 
ments (T™ 1 • • • T£ k ), with (. . .) standing for the integration 
over the joint probability density (0 and m > 0. Here, we 
calculate all the moments with J^i n i < 4 by deriving a set 
of algebraic relations for them and reducing the moments to 
forms of Selberg's integral. Presenting the relevant technical 
details at the end of the paper, we now discuss applications of 
the obtained results to various linear and nonlinear statistics 
on the transmission eigenvalues. 

Transmitted charge cumulants.- The statistics of charge q 
(in units of e) transmitted through the cavity over the obser- 
vation time is usually described by means of the current or 
charge cumulants, ((q m ))- Following Levitov et a h 11 ' 12 , it is 
convenient to use a general formula for the cumulant gener- 
ating function expressed in terms of the transmission eigen- 
values as follows £ m ^((q m )) - (£< ln[l + T t .(e x - 1)]}. 
One finds that the first cumulant, ((q)) = n(Ti), gives the 
conductance, the RMT average of which is known 4 to be 



«9» = (9) 



N X N 2 



N = Ni + N 2 



(4) 



while ({q 2 }} = n((Ti) - (Tj 2 )) yields shot-noise as follows 6 



((g 2 ))^( J ,} = ( 5 >. 



(AT 1 -l + |)(iV 2 -l + |) 



(iV-2+f)(iV-l+4) 



(5) 



An equivalent to RMT derivation of these and related results 
within a semiclassical approach may be found in Ref. 13. 

The charge cumulants are an example of linear statistics on 
the Ti's that is fully determined by the transmission eigen- 
value density, p{T) = (][]. 5(T - T,)}. However, the latter 



2 



is analytically known only in some limiting cases of a few 414 
or many 15 open channels, restricting the use of p(T) for the 
calculation of ((q m )) and full counting statistics 16,17 to these 
cases. 

In the general situation of arbitrary N12, one can alterna- 
tively consider the joint probability distribution (fJI and ex- 
ploit its simple algebraic structure (i.e. that of the Selberg's 
integral kernel) to derive exact relations for its moments. 6 For 
example, Eqs. ( TTBT l and ( TT6I ) presented below yield straight- 
forwardly and in a uniform way exact results ©, © and (0 
for (<?}, (p) and var(<7), respectively. This approach was re- 
cently extended further to find the third cumulant {{q 3 )) = 
((q— (q)) 3 } exactly. 8 For completeness and later use, we write 
down this result (following from Eq. ( TT8l below) as follows 



vax(q) (N — 3 + 
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As to the fourth cumulant of the transmitted charge, its ex- 
plicit expression can be found from Eqs. ( fT5l ) - ( fT9l ) according 
to ((q 4 )) = n[(T 1 )-7(T 1 2 )+12(T 1 3 )-6(T 1 4 )], the result being 
too lengthy to be reported here. In the case of the single-mode 
leads, Ni = N 2 = 1, one gets ((q 4 )) = - T § g , -i, -i at 
the values of (3 = 1, 2, 4, respectively. In the opposite semi- 
classical limit of many channels, Nip 3> 1, we represent the 
outcome of our calculation as the following -k expansion: 

((q 4 )) Nf - 87ViiV| + l2N?Ni - 8N?N 2 + N$ 



var(q) 



N 4 
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The leading order term agrees with the result obtained by a 
different method in Ref^£ The next order term gives a week 
localization correction which vanishes at p = 2 or Ni = N 2 . 
In the latter case of symmetric cavities, one further finds 
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Conductance cumulants.- As is clear from the discussion, 
the presented method is equally applied to nonlinear statis- 
tics determined by different transmission eigenvalues as well. 
The simplest example of such a quantity is the variance of the 
conductance, the exact RMT result of which readsi 6 - 



2(JVk 



1 



var(g) ^ 

(g) p( N -2 + Z)(N 



1 



21 

8 1 



(9) 



We note that var(g) = 2(g){p) / PN\N 2 makes a relation of 
(0 to the linear statics (01 and (f5]l considered above. 

The third cumulant of the conductance, the so-called skew- 
ness, can be also found from Eqs. (TT~5T > - ( fT8l and be repre- 
sented after some algebra in the following compact form: 
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One finds immediately that ((g 3 )} = 8[(g) / /3A?iiV 2 ] 2 ({g 3 }). 
It is also worth noting that the skewness vanishes for symmet- 
ric cavities (Ni — N 2 ) at P = 2. This holds generally for any 
odd cumulant of g (or q), as the corresponding distribution 
becomes symmetric around § (or i) in this caseJ£ 

By our method, we have also calculated the fourth cumu- 
lant ((g 4 )) which is related to the conductance kurtosis. The 
fourth moment of the conductance is found to be determined 
by moments of T„'s as (g 4 ) = n{(T?) + (n - 1)(3(T 1 2 T|) + 
4<TiT|» + 6(n - l)(n - 2){T l T 2 T 2 ) + (n- l)(n - 2)(n - 
3)(TiT2T3T4)], the corresponding cumulant being given by 
the standard formula. Since the resulting explicit expression 
appears to be too cumbersome, we restrict our consideration 
to the limiting cases of the single-mode and many-mode leads. 
In the former case, we get {{g 4 )) = -Jfg, - T | 5 , -gjg 
at P = 1, 2, 4, respectively, whereas in the latter case of 
-^1,2 > lwe arrive at the following expression: 
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One can readily see that higher cumulants contribute in the 
next order of -jj, thus the conductance distribution gets more 
Gaussian-like as N grows. 19 This tendency becomes even 
stronger for symmetric cavities at P = 2, as then ((<? 3 )) = 
identically and ((<? 4 )) vanishes in the leading and next-to- 
leading orders. In this case of N12 = n>l, one gets 
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Shot-noise variance- Now we discuss statistics of the shot- 
noise power. Its average value is given by (0), whereas its 
second cumulant, the variance, is determined by var(p) = 
n[(T 2 ) - 2{T 3 ) + <T 4 )] + n(n ~ l)[(TiT 2 ) - 2(T 1 T 2 ) + 
( T i T 2)} ~ n K T i) ~ ( T i )] 2 - However, as in all other cases 
of the fourth order cumulants considered above, the explicit 
result for var(p) can not be cast in a compact form, so we 
present the limiting cases again. In the single-mode case, we 
get var(p) = j^, correspondingly at p = 1, 2, 4. In 
the many channel case, we obtain the following expansion: 



var(p) 
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(13) 



As usual, the first weak localization correction vanishes for 
unitary symmetry, P — 2. The next order correction can be 
also found and reads as follows 
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Mehta's book 3 which read as follows 



n, 



(T 1 T 2 ---T m ) 



n 



a+|(n-i) 



l + f(2n-i-l) 
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[a + 1 + g(n - l)]ni - f (n - l)n 2 
a + 2 + /3(n-l) 



(15) 
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To calculate moments containing higher powers of T m , one 
may note that 



{T X T 2 . . . T£) = (T X T 2 . 



Q T k+1 
_ m 

9T m fc + 1 



and employ a partial integration here (for a > 0). This yields 
a contribution (TxT 2 . . -)'n-i at trie u PP er boundary T m = 1, 
where notation (. . has been introduced for an averaging 
over a joint density of eigenvalues Ti, . . . , T n _i which con- 
tains in addition a factor n"^ 1 |1 ~ This case is also 
contained in the general form of Selberg's integral. In partic- 
ular, the corresponding analogue of ( fT5l ) is found to be 



!C = (T 1 ...T m >;_ 1 = JJ 

3=1 



a + §(n-j- 1) 
+ l + f(2n-i-l) 



(17) 



FIG. 1 : The variance of the shot-noise power in chaotic cavities as a 
function of the channel numbers in the leads: a) The case of symmet- 
ric cavities (N\ — N2 = n) for three RMT ensembles, when var(p) 
saturates at the universal value (64/3) _1 at large n. b) The case of 
asymmetric cavities with preserved time-reversal symmetry (/3 = 1) 
for fixed number Ni of channels in the one lead and varied one N2 
in the other lead. The shot-noise variance shows a well-developed 
maximum at N2 ~ Ni and then diminishes down to zero according 
to var(p) » 2Ni (N1-I + j)/PN% as N 2 grows. 



In this way we are able to calculate all the moments up to the 
fourth order. For moments of the third order we have 



a+f(n-l)(l-2(T 1 lf» 



a + 3 + /3(n-l) 



(TiTi) = 



[a + /3(n-l)]n' 1 -f(n-2)n. 



a + 2 + j3(n-l) 
The fourth order moments are given as follows 



(18a) 
(18b) 



where we have omitted the terms ~ n~ 3 and put N x = N 2 for 
simplicity. The general dependence of the shot-noise variance 
on the channel numbers in the leads is illustrated on Fig. 1 . 

Covariance of g and p.- It is also instructive to consider 
statistical correlations between the conductance and the shot- 
noise power that can be characterized by their covariance 
cov(g,p) = (gp) — (g)(p)- This quantity involves moments 
of Tj's up to the third order. Surprisingly, the resulting exact 
formula for cav{g ) p)/vax(g) turns out to be given precisely 
by the r.h.s. of Eq. ©. This expression vanishes in symmet- 
ric cavities at f3 = 2 that can be easily seen again as a con- 
sequence of the symmetry of the distributions. In this case, 
therefore, g and p become uncorrected on the level of their 
averages and it would be interesting to understand to which 
extent this holds for their higher moments in general. 



Moments of {Ti} and Selberg's integral - Finally, we dis- 
cuss the derivation of general moments (T" 1 • • ■ T? k } at arbi- 
trary positive a and (3. Moments with all m — 1 as well as 
(Ti) can be found from recursion relations already given in 
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(19a) 



(19b) 



(19c) 
(19d) 

(20) 



At last, expressions (fT~5T > and (Tl\ taken at m — 1, . . . , 4 make 
the above algebraic system of equations be closed. 
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Discussion - We have applied essentially the theory of Sel- 
berg's integral to problems of quantum transport in chaotic 
cavities. The cumulants up to the fourth order for current 
and conductance fluctuations and up to the second order for 
shot-noise have been calculated exactly at arbitrary channel 
numbers and symmetry parameter j3. We note that the pro- 
posed method can be also used for determining the corre- 
sponding distributions in closed form suitable for an analytic 
work in the case of a few open channels as well as for nu- 
merical implementations. 7,20 Our analysis of shot-noise statis- 
tics suggests that in close analogy with universal conductance 
fluctuations, 19 we may characterize universal fluctuations of 
shot-noise by their cumulants as follows 

((p m » ac (p) 2 ~ m . (21) 

In the limit of large (and equal) channel numbers, Eq. d2TT > 
yields the Gaussian distribution which is peaked at (p) i=s j 
and has a width given by the universal value of var(p) » — j-^ . 
It would be highly interesting to understand how our findings, 
which are relevant to the zero-dimensional (RMT) case, can 
be extended to higher dimensions where some analytical re- 
sults are already known. 2 =i^ 

It would be also desirable to compare our results with the 



relevant experimental data. For example, measurements of 
up to the fifth cumulants of the transmitted charge have been 
recently done in a weakly coupled quantum dot. 23 However, 
performing similar experiments for a lower impedance device, 
such as a perfectly open chaotic cavity with several modes in 
the leads, requires a much higher detector resolution that rep- 
resents a current experimental challenge. On the other hand, 
the conductance distribution in chaotic cavities 24 has been 
also tested by means of experiments with microwaves^ A 
direct comparison with our results in this case could be, how- 
ever, not so straightforward, as it involves taking into account 
effects of dephasing and absorption^i 2 ^ as well. 

As another potential application of our results, we men- 
tion quantized transport in graphene p-n junctions that 
has been very recently studied both experimentally 2 - 7 - and 
theoretically. 28 The exact RMT results for higher cumulants 
of the conductance and noise may be useful for understand- 
ing possible mechanisms of edge mode mixing in the bipolar 
regime leading to chaotic transport there. 28 Further work in 
this and in the other directions mentioned above is needed. 

We thank K. Ensslin for correspondence and drawing our 
attention to Ref. 23. The financial support by SFB/TR 12 der 
DFG and BRIEF grant is acknowledged with thanks. 
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